Spin injection and magnetoresistance 
in ferromagnet/superconductor/ferromagnet tunnel junctions 
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We theoretically study the spin-dependent transport in a ferromagnet/superconductor/ferro- 
magnet double barrier tunnel junction. The spin-polarized tunneling currents give rise to spin 
imbalance in the superconductor. The resulting nonequilibrium spin density suppresses the 
superconductivity with increase of the tunneling currents. We focus on the effect of asym- 
metry in the double tunnel junction, where the barrier height of the tunnel junction and the 
spin-polarization of the ferromagnets are different, on spin injection, and discuss how the super- 
conductivity is suppressed in the asymmetric junction. Our results explain recent experimental 
results on the critical current suppression in high- Tc SCs by spin injection. 



I. INTRODUCTION 

Spin-polarized tunneling plays an important role in the 
spin-dependent transport of magnetic nanostructures M . 
First the spin-polarized tunneling causes a large magne- 
toresistance in ferromagnetic single tunnel junctions [g ; 
the tunnel resistance decreases when the ferromagnetic 
moments are aligned in a magnetic field Q. Second the 
spin-polarized tunneling current driven from ferromag- 
nets (FM) into a normal metal (N) or a superconductor 
(SC) creates a nonequilibrium spin polarization in N or 
SC |j,y. Recent experiments have shown that a strong 
suppression of superconductivity occurs by injection of 
spin-polarized electrons in tunnel junctions consisting of 
a high- Tc superconductor and a ferromagnetic manganite 

A double tunnel junction containing SC sandwiched 
between two FMs (FM/SC/FM) is a unique system to 
investigate the nonequilibrium phenomena of spin and 
charge imbalance in SC caused by the tunneling currents 
but also the competition between superconductivity and 
magnetism induced by spin polarization in SC. In a sym- 
metric double junction, where the tunnel barriers and 
the ferromagnets are the same, we have predicted an in- 
triguing magnetoresistive effect; in the antiferromagnetic 
(AF) alignment of magnetizations, the spin density accu- 
mulated in SC strongly reduces the superconducting gap 
A with increase of tunneling currents, while in the ferro- 
magnetic (F) alignment there is no such effect because of 
the absence of spin population in SC |11[| . In this paper, 
we take into account the asymmetry in the junction, and 
discuss how the spin density is accumulated in SC and 
suppress the superconductivity of SC, depending on the 
difference in the tunnel resistance of the barriers and in 
the spin polarization of FMs. 

II. FORMULATION 

We consider a FM1/SC/FM2 double tunnel junction 
as shown in Fig. 1. The left and right electrodes are made 



of different ferromagnets and the central one is a super- 
conductor with thickness d. The magnetization of FMl is 
chosen to point up and that of FM2 is either up or down. 
In the asymmetric tunnel junction, the height of the tun- 
nel barriers and/or the strength of the ferromagnets are 
different, which are characterized by the different values 
of the tunnel resistance and those of the spin-polarization 
in the junction. 

We calculate the tunneling current using a phenomeno- 
logical tunneling Hamiltonian. If SC is in the supercon- 
ducting state, it is convenient to rewrite the electron op- 
erators flkcr in SC in terms of the quasiparticle operators 
7kcr using the Bogoliubov transformation 



flkT = "k7kT +VklLki, 



where lukl 
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(1 + Ck/-£'k) with the quasi- 
A2 of SC, 



particle dispersion i?k = vCk + ^^ °f ^^' ^k being the 
one-electron energy relative to the chemical potential 
which is chosen to be zero and A being the gap param- 
eter. Then, using the golden rule formula, we obtain the 
spin-dependent currents /jo- across the ith junction: 
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FIG. 1. Double barrier tunnel junction consisting of two 
ferromagnets (FMl and FM2) and a superconductor (SC) sep- 
arated by thin insulating barriers. 



Here, Gi„ [i — 1,2) is the tunnel conductance of the ith 
junction for electrons with spin a if SC is in the normal 
state, and is given by d^ oc \T^\'^VsV'^^, where \T^\'^ is 
the tunneling probability of the zth junction and I?s and 
Ppj are the spin-subband densities of states in SC and 
FMi, respectively. The quantity Hi is given by | [l^ 



M = 



5 EN 



E^-eV^)-foiE^ 



eV,] 



(2) 



where /q is the Fermi distribution function of thermal 
equilibrium in FM and Vi the voltage drop at the ith 
junction (Vi + V2 = V). The quantities S and Q* are 
quasiparticle spin and charge densities in SC and are de- 
fined by 



5= -^(/kT -/ki), Q* ^^q-kha-, 



(3) 



ko- 



where /kcr = (7k(j7kcr) is the distribution function of 
quasiparticles with energy E\^ and spin a and gk = 
|wkP — It'kP is the effective charge of a quasiparticle in 
the state k. 

The conservation of total charge Qtot = SkCT('^kCT'^kcr) 
in SC gives Ii^+hi = h'i+hi, which yields the relation 



giPi + 92P2 S + Q*/2 = giNi - 92M2, 



(4) 



where gi — Gi/{Gi + G2) (51 + 52 = 1) is the reduced 
conductance of ith junction, and 

Pi = (Gi| — Gil) /^ii -P2 — (G21 — G2i) /G2, 

where P2 — P2 for the F alignment and P2 = — P2 for 
the AF alignment of magnetizations. Pi and P2 are the 
degree of spin-polarization of FMl and FM2. 

The quasiparticle spin density S generated in SC is 
calculated by balancing the spin injection rate (dS/dt)^. 
= [(/i| — /il) — (/2|— /2x)]/2e with the spin relaxation rate 
S/ts, where T5 is the spin-relaxation time. The result is 



S 



gig2{Pi - P2 



1 - (.91 A + 92P2Y + Ts 



{Afi+U2).. 



(5) 



where F^ = gig2{Tt/Ts), n = 2e'^Vs{Ri + R2) is the tun- 
neling time, and Rj = l/Gj. Note that the spin density 
in SC is proportional to the difference {Pi — P2) for the F 
alignment and the sum {Pi + P2) for the AF alignment. 
The quasiparticle charge density Q* is obtained by bal- 
ancing the injection rate {dQ* / dt).^. with the relaxation 
rate Q* /tq* |1^, where tq* is the charge relaxation time 
pjl, and using Eq. (H) in the form 



-(f)?| 



.giMk - .g2AA2k 



-(51^1 +52-P2)(/kT -/ki 



(6) 



where M,^ = (l/2)[/o(^k - eV,) -/o(^k + eV,)]. 

The superconducting gap A in SC is determined by 
/kcr through the BCS gap equation 



1 



^BCS 



E 



1 ^ /kT ~ /kj 



(7) 



It follows from Eqs. (|g) and (|g) that, if the junction is 
symmetric, both S and Q* vanish for the F alignment, 
while S ^ Q and Q* = for the AF alignment. In 
the asymmetric case, S and Q* become finite for both 
alignments. In the following, we restrict ourselves to the 
case TQ* <ti Tt <^ Ts, where the charge imbalance is very 
small {Q* ~ 0), so that the nonequilibrium effect is dom- 
inated by the spin imbalance. In addition, the thickness 
of SC d is much smaller than the spin diffusion length 
h = VDts, D being the diffusion constant, so that 
the distribution of quasiparticles is spatially uniform in 
SC. Then, the distribution function /to- is described by 
/o, but the chemical potentials of the spin-up and spin- 
down quasiparticles are shifted oppositely by 5^s from 
the equilibrium one to generate the spin density; 



hi ^ fo{Ek-S^is), hi = fo{Ek + S^ls)■ 



i8) 



We solve self-consistently Eqs. (||) - (^ with respect to A, 
Sfis , and Vi , and obtain A and S as functions of V. The 
results are used to calculate the total tunneling current 

linj = li] + hi '■ 

^ 1 / M + AA2 \ 1 - {giPJ + g2Pi) + Ts 

'"^ eVs\Ri+R2) 1- {giPi + .92^2)^ + Fg ' 

which we call the injection current. 

III. RESULTS 

We briefly discuss the tunnel magnctorcsistance 
(TMR) in the normal state (T > T^), in which Wi +7V2 = 
VseV, so that the TMR ratio, AR/Rp = {Ra-Rf)/Rf 
has the form 



AR 

'r^ 



4.9152^^2 



l-(.giFi+ 52^2)2 +F5 



(10) 



The TMR is degraded in the case of strong asymmetry 
in the conductances (Gi ^ G2 or Gi ^ G2). A large 
TMR ratio is obtained when the following conditions are 
satisfied; the tunnel barriers are similar {Ri ~ R2) and 
the spin relaxation time in SC is long compared with 
the tunneling time (rt/rg < 1). The latter condition 
is {pn/Ri) + {PN/R2) > i^d/l'g), where pN is the re- 
sistivity of SC in the normal state and A the junction 
area, which requires a low junction resistance and/or a 
thin SC with d much smaller than Is- If these condi- 
tions are satisfied, we have the optimum ratio AR/Rp ^ 
P1P2 1 {^ ~ P1P2) in the normal state. 
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FIG. 2. Dependence of the critical current Jc on the in- 
jection current /inj for the spin polarization Pi = 0.8 of FMl 
and different values of Pi of FM2. The open circles indicate 
the critical current measured at T = 80 K and T = 84 K 
(Tc ~ 89 K) in a FM/SC/N junction made of a high-Tc SC 
and a ferromagnetic manganite with P ~ 100% |7| . 



injection current at temperature T jT^ — 0.9 for P^ — 0.8 
and three values of Pi — 0, 0.4, and 0.8. Other parame- 
ters are taken to be Qi = 1/2 (i?i — R2) and Fg = 0. In 
the case that FMl and FM2 are the same ferromagnets 
(bottom panel), the critical current Jc in the AF align- 
ment steeply decreases and vanishes at a small value of 
/inj , whereas Jc in the F alignment shows no dependence 
on /inj . In the case that FMl and FM2 are different (mid- 
dle panel), the critical current decreases with increase of 
injection current in both alignments but in different way; 
Jc decreases more slowly in the F alignment than in the 
AF alignment. If one of the ferromagnets, FM2, is re- 
placed by a normal metal (N), we have a heterostructure 
junction FMl/SC/N, which corresponds to the junction 
with P2 = (top panel) . The calculated result for P2 = 
explains the critical current suppression by spin injection 
observed in the heterostructure junctions consisting of a 
high-Tc SC and a ferromagnetic manganite with P ^ 100 
[6-10]. 
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When temperature T is lowered below T^ at zero bias, 
SC becomes superconducting. As the injection cur- 
rent /inj increases, the superconductivity is strongly sup- 
pressed by the pair breaking effect due to the increase of 
the quasiparticle spin density S in SC. The amount of 5 
accumulated in SC is directly connected to the injection 
current /inj by the relation 



5 = 



P1-P2 



l-{giPi + g2Pi) + Ts 



eVs 



Gi 



<^2 



(11) 



When FMl and FM2 are the same (Pi = P2), the in- 
jected spins vanish in the F alignment, and are accu- 
mulated only in the AF alignment. Therefore the pair 
breaking effect occurs only in the AF alignment. How- 
ever, when the FMs are different, the injected spins are 
accumulated in proportion to (Pi — P2) for the F align- 
ment and (Pi + P2) for the AF alignment, and thus we 
have the pair breaking effect in both alignments. 

The suppression of the superconducting gap A by spin 
injection is detected by measuring the critical current Jc. 
According to the Ginzburg-Landau theory, Jc is propor- 
tional to A^, because Jc ex A^Uc, Wc being the critical 
superfluid velocity, and Vc cc A Q. Figure 2 shows the 
cube of the normalized gap, (A/Aq)'^, and thus the nor- 
malized critical current, (Jc/Jco), as a function of the 
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